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Introduction. Assuming a classical statistical system of point particles the funda- 
mental equations of continuum thermomechanics (continuity equation, equation of 
motion, and energy equation) shall be derived exactly. The macroscopic state func- 
tions (density, velocity, stress, energy density, heat flux) are interpreted as expected 
values. 

The systems can represent any continua (gases, liquids, solids). It is not assumed that 
the particles are identical. Also, the nature of the interacting forces can be different 
for different pairs of particles. Therefore it is not necessary that the particles are 
molecules in the sense of chemistry. The theory is also valid when the continuum is 
viewed as the system of its atoms (or even elementary particles). It appears to us 
physically justified to assume that atoms (or at least the elementary particles) are 
formed of points and rotational and internal degrees of freedom can be neglected. 
Hence, for the general theory there is no difference between a mixture and a chemical 
compound. The difference is solely due to the nature of the mutual potentials. The 
failure of classical mechanics in the atomistic realm and the influence of quantum 
mechanical effects, however, are not taken into account in our considerations. 

In section 2 we will formulate the problem precisely and summarize the results. By 
introducing appropriate abbreviations in section 1 the considerable difficulties in as- 
signment will be eased. In section 3 the sufficient conditions will be formulated under 
which the following investigations are valid. Sections 4 and 5 contain the derivation 
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of the expected values of the state functions. In sections 6 and 7 we will prove the 
validity of the fundamental equations, first for the case that no external forces are 
present. In section 8 we will discuss the influence of external forces. Our investi- 
gations are crucially based on two mathematical Lemmas, formulated and proved in 
section 9. 

The problem being treated here was tackled by Irving and Kirkwood [lj. This work 
differs from pQ in the following points: 

1) The proofs satisfy all requirements of mathematical rigor and, additionally, 
should be easier than the ones in pQ. 

2) For the stress tensor and the heat flux we give closed integral expressions. In 
[T] these quantities are presented as infinite series which only makes sense if the 
probability density as a function of the spatial variable is analytical everywhere. 

3) The interpretation of stress tensor and heat flux as expected values will be 
carried out in detail. This is only attempted in pQ. 

4) The assumption that all particles are identical is not made. 

5) The external forces do not have to have a potential and may — apart from space 
and time — depend also on the velocities of the particles. 

6) The use of 5 functions has been avoided. In pQ they only serve technical issues. 

1 Definitions and assignments 

a. Vectors and tensors. Vectors and points (point vectors) are described by lower 
case letters in bold font, tensors of higher order are described by capital letters in bold 
font. The product of two tensors (dyad) will be represented by the symbol <g>, the 
inner product of two vectors or of a vector and a tensor will be represented by a dot 

and V x F(x) represents the differential (gradient) of the tensor F(x). V x • F(x) 
represents the differential (divergence) of F(x). 

b. Probability density. We consider a system of point particles. The j th parti- 
cle is (j), its mass is mj. A state of this system is characterized by the par- 
ticle locations xi,--- , xjv and velocities - , £/v, indicated by a "point" 
(xi,--- , Xjv;£i,-" , &v) = ( x i;£i) of 6iV-dimensional phase space fL This phase 
space is the product of 2N three dimensional Euclidean vector spaced. The proba- 
bility density of the states in 0, at time t is denoted by 

W = W(xi,~ - .xjv;^,-- - ,&v;t) = W(xi;&;t). 

c. Potential energy. We assume that the force exerted by (j) on (k), kffc, depends 
only on the locations xj and x^. By invariance it follows that \s.jj~ must be a central 
force whose contribution is dependent only on the distance = \xj — x^|. Newton's 

The Xi and vary in the entire, infinite three dimensional space and not only in a subregion of 

it. 
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third law says that kj k = — kkj. The particle pair (j),(k) therefore corresponds to a 
function, 

(1-1) Vjk (r) = V kj (r), 

such that Vjk{rjk) yields the potential energy of the pair. Furthermore, 

k jk = -V Xj V jk (r jk ) = -V! k (r jk ) X3 * fc 

(1.2) jk 

= -'kkj = -V Xk V jk (r jk ) = V^ k {r jk )— K 

I jk 

For the total inner potential energy U of the system this results in 
(1-3) U = C/(xi,-- - ,xtv) = ^2v jk (r jk ) = y^V jk (r jk ). 

j<k j^k 

Here, one has to sum of all those pairs (j, k) for which j < k or j / k is true. It is, 

N N _ 

(i-4) v x r ^v x v jk {r jk ) ^r/,(/>) X " x " 

fe=l k=l 
k^j k^j 



d. External forces. Apart from those forces acting on particle (j) due to other particles 
we assume that yet another external force kj is exerted that depends only on location 
Xj and the velocity £j of the particle at time t, 

kj = kj (xj , £j , t) . 

Restrictively, we require, however, that the functions kj obey the equation 

N 1 

(1.5) ^ v £ •k ; ix,.£,./i ,) 

identically. In particular, this condition is met if the kj(xj, £j, t) do not depend on 

is- 

e. Integrals. The (6N — 3)-dimensional subspace of S7, when discarding the spatial 
variable Xj belonging to (j), shall be denoted Qj. Analogously, we denote the subspace 
Cljk that comes from discarding Xj and x^. Let F be a scalar, vector, or tensor function 
defined over 0. Then, the abbreviations, 

(F\xj = x), (F\xj =x,Xjk = y), 

mean that F is the integral over flj (0,-fc, respectively), where, after performing the 
integration, the free variable Xj (or the variables Xj and x&, respectively) are to be 
replaced by x (or by x and y), respectively. Generally, we call 
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/(y) dy 

'y 

the volume integral of f(y) over the infinite three-dimensional space of locations y. 
Then, apparently, the relation 

(1.6) / (F|x, = x, x fc = y) dy = (F| Xj = x) 

is true. 
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2 Posing the problem 



We assume that the probability density, W(xj;£j;i), is defined for all (xj,£j) E ^ 
and is continuously differentiable with respect to all variables. Under the restriction 
that the external forces fulfill the condition (|1.5|l ^l the principle of the conservation of 
probability in phase space yields the differential equation in conventional fashion, 

(2-i) ^ = E {"6 ■ + ~ k *) • V ^ W } > 

which determines the rate of change of W. 

It is the task of this study to derive the fundamental equations of thermomechanics 
from (|2,ip under the regularity requirements that will be formulated in section 3. 
First we will assume the absence of external forces so that (12. ip is simplified to 

TO ^ = E ■ + ^ v *^ • • 

The case of kj ^ will be treated in section 8. For k, = the fundamental equations 

are as follows: 

A) Continuity equation 

(2.2) | + V X - (H = 0. 
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B) Equation of motion 



(2.3) p ^ + u.V x u^ =V X -S. 

C) Equation of energy 

de 

(2.4) _ + Vx .(q-S-u + eu) =0. 

These fundamental equations connect the following macroscopic state functions: 

p = p(x, t) = Mass density, 
u = u(x, i) = Velocity, 
S = S(x, i) = Stress tensor, 
e = e(x, i) = Energy density, 
q = q(x, t) = .ffeai flux density. 

We will prove that these fundamental equations are true if the aforementioned state 
functions are interpreted as expected values. In sections 4-5 we will show that these 
expected values are given by the following expressions, 

(2.5) p = J2m j (W\ x i= x), 

3 

(2.6) pu = E m i^i^ I x i = x ) - 



2 This requirement ensures that phase space is locally volume preserving. 
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The stress tensor is symmetric and consists of a kinetic contribution, Sk, and a 
interaction contribution, Sy, 

(2.7) S =S K + S V , 

(2.8) S K = - m i<& " u ) ® & " u ) w I x i = x )> 

3 

(2.9) S v = 

I \ ~TT~ V jk(\ z \) I (W \ Xj = x + az,x k = x - (1 - a)z)da\ dz. 
1 i +k J* V l z l J<*=o ) 

The energy density splits into the kinetic energy density, ex, and the interaction 
energy density, ey, 

(2.10) e =e K + e v , 

(2.11) e K = ±Y i m j ($W\x j = x), 

3 

(2.12) e v = \ J2( V 3k(\*3 - x k \)W | xj = x). 

The heat flux has three terms, the kinetic contribution, qK, the transport contribution, 
qx, and the interaction contribution, qy, 

(2.13) q = qx + qT + qv 

(2.14) q K = \Y, m Mi ~ u | 2 (^i " ™)W | x, = x), 

j 



(2.15) q T = \ Y^Ui ~ u)^fc(|xj- - x fe |)W | Xj - = x), 

(2.16) qv = "^E / ||| z -^(l z l) 

• j / f ^ J ^ — W | Xj = x + az, Xfc = x — (1 — a)z^ da| dz. 



3 Regularity conditions 

The expectation expressions occurring in (|2.5[) - (|2.16|) are improper integrals over W 
and Vjk- It is therefore clear that certain regularity conditions must be met for W 
and Vjk if the state functions in (|2.5|) - (|2.16|) are well-defined and the functions of x 
and t are continuously differentiable. Already the condition J n Wd£l = 1 requires 
that W approaches zero sufficiently fast as |x,-| — > oo and — » do. 

The following three conditions are sufficient for the validity of the results in this work: 
A) There is a number 5 > such that the function 

N N 

Gfo&t) = Wix^t) [] |x/ +<5 [] |£ fe | 3+<5 , 

j=i k=i 
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as well as its derivatives are restricted by a constant solely dependent on t. 

B) The functions Vj k (r) are defined for all r, continuously differentiable and, together 
with their derivatives, are finite. 

C) The functions k, (x, £,t) are defined for all values of x, £,t, continuously differen- 
tiable, and with the constants A(t) and B(t), solely dependent on time, and satisfy 
have, 

|k,-| < A(t)\£\ + B(t), |V € k*| < A(jt)\£\ + B(t). 

These conditions are sufficient for the convergence of all the improper integrals, in- 
terchanging the order of integration, and for differentiation and integration etc. 

Furthermore, condition A ensures the validity of the following Lemma, which is proved 
by partial integration: 

Suppose that F(x;£i) is a continuously differentiable function defined in f2 which, 
with the constants A and B, satisfies the inequalities, 

N N N 

\F\ <AY[\£ k \ 3 + B, |V Xj F| <Al[\t k \ 3 + B, \V^F\<All\t k \ 3 + B. 

k=l k=l k=l 

Then, these formulas are valid: 

(3.1) J FV Xj W = - J WV^F, J FV^W = - J WV^F, 

where one integrates over a subspace of £1 that contains the space of Xj or £j, re- 
spectively, as another subspace (omitting the volume elements). In particular, if F is 
independent of Xj or £j, one has, 

(3.2) jFV X:j W = or J FV^W = 0, respectively. 



Note I. The equation (|2.ip is a linear partial differential equation of first order for W. 
According to generally known rules W(xj;£j;i) is therefore uniquely determined for 
all t if the initial density, W(xj;£j;0) = Wo(xjj^j) is prescribed. Therefore it would 
be desirable to apply regularity conditions only to Wo( x «;&) and not to W(xi',£i',t), 
and then prove A as a property of W. The author, however, did not yet succeed in 
this. 

Note II. If W is not differentiable for all (x$; G Q, and t the equation (|2.1[) generally 
loses its meaning. The statistical mechanics, however, is meaningful if W is only 
integrable over 0. It is possible to generalize our investigations to this case. The 

N 

condition A can then be replaced by the requirement that W \\ |£fc| 3 be integrable 

k=l 

for all t over fL Then, the expressions f)2 .5^ — f|2 . 16j) remain meaningful. The state 
functions, however, being functions of x and t are no longer continuously differentiable 
for all values of x and t such that the fundamental equations in the form (|2.2p - (|2.4p 
in general, have no meaning. They must be replaced by the laws of conservation of 
mass, momenta, and energy in finite form. 

To carry out this sort of generalization demands considerable formal work. But this 
can be avoided by understanding all the present differentiations in the sense of the 
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theory of distributions by L. Schwartz [2j. Then, (|2.ip is completely equivalent 
with the principle of conservation of probability. The fundamental equations (j2.2j) — 
(|2.4p are generally valid in the sense of the theory of distributions and also when 
the probability density W does not even exist. In this latter case W and the state 
functions are to be viewed as measures. The fundamental equations are valid in 
the conventional sense only for those values of x and t for which the state functions 
represent continuously differentiable functions. One can readily derive the transition 
conditions at the discontinuity surfaces (collision and acceleration waves) once (|2.2p ~ 
(|2.4p are understood in the sense of distributions. 

4 Expected values of the state functions 

The value of a physical quantity, F, for the particle (J) shall be given by a function, 
fj(?ti',£i)t defined in 0. According to the rules of probability calculus and according 
to section le one has 

(4.1) £j(F) dx = (fj W\x j = x)dx 

as the expected value of F for (j) under the condition that (j) is in the volume element 
<ix at x. Hence, the density £ (F) of the expected value of F for all particles at the 
position x is 

(4.2) £(F) = = Y^fJjW | Xj = x). 

3 3 

The expected value of F for (j) under the condition that (j) is in the volume element 
<ix at x and (&) in the volume element dy at y is given by 

(4.3) £jk(F) dx dy = (fj W \ xj = x, x k = y) dx dy. 

a. Density and velocity. It follows immediately from (|4.2p that (|2.5p and (|2.6p repre- 
sent the expected values for mass and momentum densities. 

b. Stress tensor. Let J be a region of three-dimensional space with a continuously 
differentiable boundary, T . We call A the exterior of J with n x being the normal 
unit vector at x £ T pointing outwards, and with dT x the associated surface element 
(see Fig. Q}@ 

As generally known, the stress tensor is characterized by the fact that for any part 
J of the considered body, the force, k, exerted by A on J can be represented in the 
form 

(4.4) k = ^S(x)-n x d^ x . 

One readily sees that according to (12. 8j) . 

(4.5) k K = J S K • n x dF x 

translator's footnote: The image was cropped from a pdf of the original paper. Note the differ- 
ences in font styles in the figure and in the text. 
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gives the expected value of the "kinetic" force that corresponds to the momentum per 
unit time transported from A to J . Therefore J has to be thought of moving along 
with the average velocity u. The kinetic contribution (|2,8p corresponds to what is 
normally called "viscous tension" , in the kinetic theory of mono atomic gases. 




Figure 1: 

Now, let us assume that the particle (j) is at the position Xj = v G J and that the 
particle (k) be at the position Xfc = w 6 A- Then, according to (|1.2j) . (k) exerts the 
following force on (j), 

v — w 

(4.6) k ifc (v, w) = -V v F jfc (|v - w|) = -VU\v - w|), , . 

v — w 

According to (|4.3|) the expected value of the force, exerted by (k) in A on (J) in J , 
is therefore given by 




= v, Xfc = w) dv dw. 



After summation over j and k we obtain from the expected value ky of the total force 
that is exerted by the particles in A on those in J . With (|4.6p this yields 

k v = -/ / v/t^^^Iv-wIKWIx^v.x^w)) dwdv. 
JvejJvreA^ l|v-w| 1 J 

As one can readily see, the integrand fulfills the conditions A, B and C of section 9. 
According to Lemma II, (19.4p and with (|2.9h it therefore follows that 

(4.7) k v = j^Svn x dJT x . 

The force k exerted by A on J is composed from the kinetic force kx ()4.5p and the 
interaction force ky (14. 7p . Because of (12 .7p we therefore have, 

k = k K + k v = J (S K + S v ) • n x dT^ = j S • n x dJT x . 
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This equation satisfies (|4.4p and so the expected value of the stress tensor is indeed 
given by (ZTMZSll . 

c. Energy density. Let us think of the kinetic energy, of particle (j) as being 
localized at the position Xj of this particle. Then, according to (|4.2p the expected 
value, eK(x), of the kinetic energy density at position x is given by (|2.1ip . 

The localization of the potential energy demands a certain arbitrariness as it is as- 
signed to particle pairs, and not — as for the kinetic energy — to individual particles. 
We assume that the potential energy corresponding to the pair (j), (k), V^(|xj — x k \), 
is distributed equally at the positions x,- and x^. Then, at position x,-, the total po- 
tential energy is, 

1 N 

(4.8) Vj = Vj(xi,--- , xjv) = Vj(xi) = 2 X/ V jk( r jk), r jk = |xj -x k \. 

k-l 

Then, according to (|4.2p . (|2.12p yields the expected value, ey(x), of the density lo- 
calized according to (|4.8[) . 

d. Heat flux (density). The heat flux, q = q(x), is characterized by the fact that for 
any part J of the body, the energy transferred per unit time from J to A, can be 
represented in the form, 

(4.9) Q = jf q(x)-n x d^ x , 
(see Fig.l). 

One is readily convinced that with (12.14p the expected value of the kinetic energy 
flowing per unit time from J to A is represented by the expression 

(4.10) Qk = J^ K ■ n x dT x . 

Therefore, qx, given by (|2.14p . is indeed the kinetic contribution of the heat flux. 
This is well known from the kinetic theory of mono-atomic gases. 

Let us consider the expression, 

(4.11) qx = q T + evu = - V jk {\xj - x k \)W \ x, = x), 

and compare to (|2,12p and (I2.15p . According to (|4.8f) we have 

(4.12) <& = Y t (t j VjW\K l =iL). 

j 

The expression £jVj(x.j) indicates the flow of potential energy, localized according to 
P~8"j) at Xj. Therefore, with (|4.12[) and because of (|4.2p . we have 

Qt = J Qt • n x dT x 

as the potential energy that must be transferred per time unit through T . However, 
J has to be seen as fixed in time. In order to obtain the potential energy Qt 
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that is transported per time unit via the moving surface T one has to subtract the 
macroscopic convection contribution, 

Qt = J ey u • n x dF x . 

So, because of (|4.11|) one finds 

(4.13) Q T = Q*-Q° = J^q T -n x d^ x . 

Therefore, the qT given in (|2.15p is indeed the contribution of the heat flux that stems 
from the transport of the potential energy. 

Finally, let us consider 

qv = qv - s v • u 

(4-14) =-^E/{^T z -^(N) 

j \ ~ ''^2 ^ ^ I x i = x + az ' x fc = x ~ (1 ~ Oi)z^ da | dz, 
(compare to (|2.9p and (| 2 . 1 6 p V In section 5 we will show that 

(4.15) Qv = ^qy • n x AF X 

represents the expected value of the potential energy, transferred per unit time from 
J to A, which stems from particles in A performing work on particles in J . Here, 
J can be viewed as fixed in time. In order to obtain the interaction energy Qy, 
transferred from A to the moving J , one has to subtract from Qy the macroscopic 
interaction work performed by J on A per unit time, 

Q y = - J (S v • u) • n x dF x . 

So, according to (|4. 14[) one finds that 

(4.16) Qy = Q*y-Qy = J^dy + S v • u) • n x dT x = J q v • n x dT x . 

Consequently, (|2.16p indeed gives the interaction contribution to the heat flux. 

Q = Qk + Qt + Qv is the total energy per unit time transferred from the moving J to 
A, and (021 is satisfied because of fljJDJ) , (|H5|) . ([435]) . and (1243]) . It has therefore 
been shown that the expected value of the heat flux is indeed given by (|2.13|) - (j2.16p . 



5 Interaction contribution to the heat flux 

We wish to prove the validity of equation (|4.15p . replacing expression (|4.14p by q v . 
To this end we investigate the mechanical system S composed solely from particles 
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in J . For this system, the forces exerted by the particles in A on the particles in J 
are considered to be external forces. These forces have the time dependent variable 
potential, 

(5.1) V(x 4 ej;t)=^2 E Ml x i - X *(*)D- 

Here, Xj G J is to be viewed as the independent variable while Xf~{t) are those 
functions that describe the trajectories of the particles G A. The internal energy 
of the system S is Ej. Then, the total energy E with respect to the external potential 
(|5,ip is given by 



(5.2) E = E I +^(x i G J;t). 

According to the energy law in particle mechanics, 

E = J^te ej;t) 



x 3 ejx k eA 1 3 feWI 



along with (|5.1|) 

Xfc(t) 

and because of x^ = £fc 



(5-3) E = - E V jk( r jk) Xj Xfc • €*, r jk = jxj - x fe |. 

The energy, however, does not equal that energy that is localized in J according 
to (|4.8p . The energy E* localized in J according to that statement is composed from 
the internal energy Ej and half of the potential energy that corresponds to all those 
particle pairs (j), (k) for which xj G J and x k G A. Therefore, we have: 

(5-4) E*=Ei + lj2 E V 3kirjk). 

From this, (|5,ip and (|5.2p results 

E * = E ~ \ E E V Mrjk). 

Xj&J X k &A 

Differentiation with respect to t yields 

(5-5) E* = E - \ £ £ lft(r iJfc ) • & - 

x 3 e.7x fc e^ jfc 

where Xj and x k have been replaced by £j and respectively. Insertion of (|5.3p into 
(|5.5p results in 

Xj - x k £j + £ fc 



Xi€jx k £A JK 
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This expression indicates how much the energy, localized by statement (|4.8p in J', 
varies per unit time and due to the interaction between the particles. The energy flux 
contribution directed from J to A is therefore given by — E* . According to (|4.3p . the 
expected value of this contribution has the form, 




The integrand fulfills the conditions A,B,C in section 9. With (|4.14p this results 
according to Lemma II, (|9.4p in 

Qv = j qv • n x dfx- 

Q.E.D. 



6 Continuity equation and equation of motion 

The proof of the continuity equation (|2.2p is very easy. Multiply (2.1a) with mj, 
integrate over Qj, replace Xj by x and sum over j. Because of (13. 2p various terms 
cancel. Eventually, one obtains equation (|2.2p with (|2.5p and (|2.6p . 

We multiply (2.1a) with mj£j, integrate over Qj, replace the free variable xj by x and 
finally sum over j, then, taking (|3.2p and (|2.6p into account, we obtain the equation 

d 

(6-1) — (pu)=vi+v 2 , 

where 

(6.2) v x = - J2 ™i itj & • V x , W) | Xj = x) , 

3 

(6.3) v 2 = ^((V Xj [/ • V^W)^ | x, = x). 

j 

Exchanging integration and differentiation with respect to Xj in (|6.2p yields, 
vi = V x • j - X) ® ^)W I xj = x) j . 

Since, 

(£j ® £i) = (£j - u) ® (£/ - u) + u ® £j + (8 u - u ® u, 
and taking (12. 5p . (12. 6p . and (12. 8|) into account, it follows that 

(6.4) vi = V x • {Sk — u ® pu} = V x • Sk - uV x • (pu) - pu • V x u. 
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Inserting (|3.ip in (|6.3p yields 

V2 = -^(V Xj W|x i= x). 



If we replace V x ^ U by expression (II. 4ft we find 




One can easily see that the integrand fulfills the conditions A, B,C in section 9. In 
particular, the validity of (|9.ip is found if the summation indices j and k are exchanged 
in (|6.5p . According to Lemma I, (j9.2[) follows with (12. 9p . 

v 2 = V x • S v . 

We insert this and (|6.4p into (|6.ip . and obtain, 

d dp du 

— (pu) = u— + p— = -uV x ■ (pu) - pu • V x u + V x • (S K + S v ). 

Bearing in mind the continuity equation (12, 2\i the equation of motion (12. 3p indeed 
follows. 



7 Energy equation 

We multiply (2.1a) with (l/2mj£j + Vj?) (see (|4.8p . integrate over Qj, subsequently 
replace xj by x and sum over j. Because of (|3.2p various terms cancel. Eventually, 
with (j2~T0l) - (j2TT2|) we obtain 

t» -, \ d e 

(7.1) -^ = qi+Q2 + <?3, 

where 

(7-2) qi = ™J>j(£fe ■ Vx,^ I x, = x), 

j 

(7-3) q 2 = - Y.^ 1 ■ I ^ = X >' 

j I 

(7-4) q 3 = -\ YffixjU ■ V^W | x, = x). 

3 
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Exchanging integration and differentiation with respect to x,- in (|7.2p yields 
Qi = -V x • {\ m M*i W I x i = x »- 

3 

Because 

%Zi = & " u ) 2 (^i " u ) + 2 [fe " u ) ® (£j - u)] ■ u + u*j + u 2 fo - u), 
and taking (|23]l . ([21)]) . ([2^]) . ([2TT]) and ([27Li]) into account, it follows that 
(7.5) gi = -V x • (<?k - u • S K + ue K ). 

Inserting (|4.8h into (|T.3[> and bearing (]3.2p in mind, we obtain 

(7-6) g 2 = ^{ l)/,iV,/,iC, • V x , W | x, = x) 

+ (^fc(^fe)Cfe • V Xfe T^ | Xj = x)} 

According to the product rule, 

V jk (r jk )V Xj W = V Xj [V jk (r jk )W] - V Xj V jk (r jk )W. 

If one inserts this into the first term on the right hand side of (|7.6|) and rearranges 
the second term according to equation (|3.ip this yields 



(7-7) q 2 = ~\ J2{-^jV jk (r jk )W | x, = x) 

+ (£j • V Xj V jk {r jk )W | Xj = x) + • V Xk V jk (r jk )W | Xj = x)}. 
Equation (HO) , applied to (j7T4l) . yields 

(7.8) © = - • V Xj EW | x, = x) = - • V x ,F ifc (r ife )W | x, = x), 

3 j^k 

where (ll.4p has been employed. Via addition of (17. 7p and (|7.8p and by insertion of 
(11.21) we obtain 



(7.9) 52 + 53 = -V x - qx - 5o, 

where q^ is given by (|4.1ip and where 



<?0 

> 

3+k 



Because of (jl.6p , we can write 50 in the form of 
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In the same way as for (|6.5j) one easily sees also here that the conditions A, B, C 
in section 9 are fulfilled for this integrand. Therefore it follows from Lemma I, (I9.2p 
that 

qo = -V x • q v 

is true where qy is given by (|4.14p . Inserting this into (|T.9[) it follows with (j7.5|) . 
(TTll . gH]), and (I4T41) that 

de 

-Q t = - V x • {q K - u • S K + ue K + qx + Qv} 

= -V x • {qx + qT + qv - u • (S K + S v ) + u(e K + e v )}, 
i.e. the energy equation (|2,4I) , 



8 External forces 

When the external forces hj are different from zero and provided that (jl.5p is fulfilled, 
the following situation emerges: 

a. The Continuity Equation ()2.2[) remains valid. 

b. According to 

(8.1) f = f(x,t) = ^2(kjW |xj =x) 

the Equation of Motion (|2.3|) transforms into 

(8.2) p f ^ + u • V x u^) = V x • S + f 

Here, f(x, t) represents the expected value of the external force density, according to 
T2D. 



c. With 

(8.3) A = A(x, t) = ' k J W \ x i = x ) 



the Energy Equation (|2.4p becomes 

(8.4) — + V x • (q - S • u + eu) = A. 

Here, and according to (14. 2)1 . A(x, t) corresponds to the expected value of the work 
performed by the external forces per unit volume and unit time. 

To prove a we bear in mind that the equation (|2.ip is obtained by adding the term, 

N 

mi 



- 1 



i=i 
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to the right hand side of (2.1a). According to the explanations in section 6 one has 
to add to the right hand side of (2.2) the term 

-EE^w^ix^x). 

3 I 

Rearranging with the help of (13. ip . we obtain 

This, however, entirely disappears according to (|1.5|) . whereby we have proved claim 
a. Completely analogous calculations prove claims b and c. 

From (|8.1|) and (|8.3p one can see that the functional form f(x, t) and A(x,t) do 
not only depend on the the external forces kj but also on W, i.e. the respective 
microscopic state of the system. In certain cases knowledge about certain macroscopic 
averages over W suffice for the determination of f and A. We will treat two of such 
cases: A) Electrical and gravitational fields. If the external forces stem from an 
electric field of strength e(x,i), we have 

(8.5) kj(x,f,t) = e,e(x,i), 

where ej is the charge of particle (j). For the expected value of the charge density A 
and the electric current density i, we obtain according to (|4.2p 

(8.6) A = A(x, t)=^2 ej (W \ Xj = x), 

j 

(8.7) i = i(x, t) = ejitjW | Xj = x). 

j 

With (|H3ME21>, (HO) and yield: 

(8.8) f = Ae, A = i • e. 

For the case of a gravitational field, g(x, t), we have to set ej = rrij and obtain, 
because of (|2.5p and (|2.6p . 

(8.9) f = pe, A = pu-g. 
B) Magnetic field. In this case 

(8.10) \c j (x,t,t)=e j £xb(x,t), 
where b is the magnetic field strength. Because of 

V € -k i (x,€,t) = e i (rot^)'b = 

condition (jl.5p is fulfilled. For a restricted/limited b(x, t), the kj satisfies the reg- 
ularity condition C in section 3. Inserting (I8.10p in (|8.ip and (j8.3fl . and with (I8.7p 
borne in mind, we obtain 

(8.11) f = bxi, A = 0. 
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9 Two Lemmas 



Let /(v, w) be a scalar vector- or tensor function of the two vector variables v and 
w that satisfy the following conditions: 

A) /(v, w) is defined for all v and w and continuously differentiable. 

B) There exists a 5 > such that the function, 

g(v,w) =/(v,w)|v| 3 +V| 3+5 , 

as well as its derivatives are bounded. 

C) The functional relation 

(9.1) /(v,w) = -/(w,v) 
is true. 

Under these circumstances the following two Lemmas are valid. 
Lemma I. 

(9.2) //(x,y)rfy = ~V x - /{z® f /(x + az,x-(l-a)z)da}dz 

Jy Z Jz Ja=0 

Proof: The conditions A and B ensure that the absolute convergence of the improper 
integrals, occurring in the following, as well as the validity of the anticipated exchanges 
of integration sequences etc. According to (|9.ip we have, 



/ /(x,y)dy = - / /(y,x)dy. 



If we introduce the new integration variable z = x — y into the left, and z = y — x 
into the right hand integral, we find 

/ /(x,y)dy= / /(x,x-z)(fz = - I /(x + z,x)dz 

Jy Jz Jz 

(9.3) =i^[/(x,x-z)-/(x + z,x)]dz. 

According to the chain rule this gives 

V x /(x + az, x - (1 - a)z) = V v / + V w / 

and 

^/(x + az,x - (1 - a)z) = z • (V v / + V w /) 

where we have to insert on the right hand sides v = x + az and w = x— (1 — a)zas 
arguments of V v / and V w /. Therefore, we have 

z • V x /(x + az, x — (1 — a)z) = —fhc + az, x — (1 — a)z). 

da 

Integration of this equation with respect to a from a = to a = 1 yields, 
z • V x / /(x + az, x — (1 — a)z)da = /(x + z, x) — /(x, x — z). 

Ja=0 
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Insertion into yields IjOjI . Q. E. D. 

Lemma II. Let J be any region in space with piece-wise smooth border surface T . 
Let A be the exterior of J and n x the normal unit vector at point x on T pointing 
outwards (see Fig. 1). Then, 



(9.4) [ I /(v,w)dwdv 



/(x + az, x — (1 — a)z)(z • n x ) da dz dJ 7 ^. 

i F Jx Ja=0 




Proof: First, one sees immediately that because of the antisymmetry of /(v, w) (|9.1 



/ / '<*■ 

./we/7 



, w) dv dw = 0. 

iveJ Jwej 
Therefore, we have 

(9.5) j f /(v,w)dvdw = / / /(v,w)dwdv. 
JveJ JweA JveJ J-w 

Now, according to Lemma I, 

(9.6) / /(v,w)dw = V v -g(v), 

J w 

where 

(9.7) g(v) = ~J { z(g) / /(v + az,v - (1 - a)z daj dz. 
According to Gauss' theorem, 

/ V v • g(v) dv = g(x) • n x dJT x . 
From this and from (|9.5p - (|9.7p results in the relation (|9.4p claimed. Q. E. D. 
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Appendix 

The following changes were made to equations in the original manuscript. 
1 . Equation (jl.4p : replaced Vxj with V x • . 
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2. Equation ([2*13]) : replaced (£.,• - u)((£j - u) 2 with - u| 2 (^- - u). 

AT AT 

3. Section 3: replaced G(x;;&;t) = Wfax^t) JJ |xj| 3+<5 |^| 6+<5 with 



j=l k=l 

N N 

3+8 



G(xi; £; t) = Wfc; Xi; i) JJ |x, | 3+5 [] |& 

i=i fc=i 

4. Equation (pT2> : replaced FV Xj W = with / FV Xj W = 

5. Equation f|T.3[) : replaced V X2 with V Xj . 

6. Section 8: replaced ^^rrij ^ — ' | x^ = x^ with 

7. Proof of Lemma II: replaced / / /(v, w) civ <iw = with 

/ / /( v ) w ) d\r dw = 
Jvej Jviej 

8. Equation ([93]) 2: replaced / / /(v, w) ciw dv with 

Jvej JweA 

/ / /(v, w) ciw dv 
Jvej J-w 
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